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ABSTRACT 

Many galaxy clusters host megaparsec-scale radio halos, generated by ultrarelativistic electrons in the magnetized 
intracluster medium. Correlations between the synchrotron power of radio halos and the thermal properties of 
the hosting clusters were established in the last decade, including the connection between the presence of a halo 
and cluster mergers. The X-ray luminosity and redshift-limited Extended GMRT Radio Halo Survey provides a 
rich and unique dataset for statistical studies of the halos. We uniformly analyze the radio and X-ray data for the 
GMRT cluster sample, and use the new Planck Sunyaev-Zel’dovich (SZ) catalog to revisit the correlations between 
the power of radio halos and the thermal properties of galaxy clusters. We find that the radio power at 1 .4 GHz 
scales with the cluster X-ray (0. 1-2.4 keV) luminosity computed within R 500 as P 1.4 ~ L 2 ^ 0 ' 2 . Our bigger and 
more homogenous sample confirms that the X-ray luminous (T 500 > 5 x 10 44 erg s _1 ) clusters branch into two 
populations — radio halos lie on the correlation, while clusters without radio halos have their radio upper limits 
well below that correlation. This bimodality remains if we excise cool cores from the X-ray luminosities. We also 
find that Pi 4 scales with the cluster integrated SZ signal within R 500 , measured by Planck, as P 1.4 ~ Uoo > 311 
line with previous findings. However, contrary to previous studies that were limited by incompleteness and small 
sample size, we find that “SZ-luminous” Uqo > 6 x 1 0 5 Mpc 2 clusters show a bimodal behavior for the presence 
of radio halos, similar to that in the radio-X-ray diagram. Bimodality of both correlations can be traced to clusters 
dynamics, with radio halos found exclusively in merging clusters. These results confirm the key role of mergers for 
the origin of giant radio halos, suggesting that they trigger the relativistic particle acceleration. 
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1. INTRODUCTION 

The presence of non-thermal components (relativistic par- 
ticles and magnetic fields) mixed with the thermal intraclus- 

ter medium (ICM) has been revealed by radio observations of 

galaxy clusters showing diffuse, giant Mpc-scale synchrotron 
radio halos (RHs) and radio relics in a substantial fraction of 

massive clusters (e.g., Ferrari et al. 2008; Cassano 2009; Feretti 

et al. 2012 for reviews). 

Giant RHs are the most spectacular and best studied cluster- 

scale non-thermal sources. Their origin is still poorly under- 

stood. One possibility is synchrotron emission from secondary 
electrons generated by the collisions between cosmic ray pro- 
tons and thermal protons (so-called “secondary models,” e.g., 
Dennison 1980). However, the same collisions should produce 
gamma rays through the generation and decay of neutral pions. 
The non-detection of nearby galaxy clusters in the y-ray band 
at 0.1-100 GeV puts serious limits on the contribution of sec- 
ondary electrons to the RH emission (Ackermann et al. 2010; 
Jeltema & Profumo 2011; Brunetti et al. 2012). 

A second possibility is that the turbulence, generated in the 
ICM during cluster-cluster mergers, re-accelerates preexisting 
GeV electrons (e.g., Brunetti et al. 2001; Petrosian 2001). The 


“historical” motivation for turbulent acceleration for the origin 
of RHs comes from the high-frequency steepening of the Coma 
halo spectrum, implying that the mechanism responsible for 
the acceleration of the emitting electrons is inefficient (e.g., 
Schlickeiser et al. 1987). More recently, the discovery of RHs 
with extremely steep spectra, 10 a ~ 1.5-2, support turbulent 
re-acceleration and disfavor a “secondary” origin of giant RHs 
(e.g., Brunetti et al. 2008; Dallacasa et al. 2009; Giovannini 
et al. 2009; Macario et al. 2010, 2011; Giacintucci et al. 2011, 
2013; Bonafede et al. 2012; Venturi et al. 2013). 

Studies of statistical properties of giant RHs in clusters and 
their connection with the cluster dynamics are extremely useful 
to constrain the origin of halos. From the inspection of the 
NRAO VLA Sky Survey (NVSS; Condon et al. 1998) fields 
containing X-ray selected galaxy clusters, Giovannini et al. 
(1999) concluded that RHs are rare at low X-ray luminosities 
{Lx < 10 45 hi} erg s _1 ), while only the most X-ray luminous 
systems host RHs, with a probability of ~l/3. Since then, a 
number of correlations have been found between thermal and 
non-thermal cluster properties, suggesting a tight connection 
between them. In particular, the synchrotron monochromatic 


10 Here, we adopt the convention f v oc v “ . 
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radio power of halos at 1.4 GHz (Pi. 4 ) has been found to 
increase with the cluster X-ray luminosity, temperature, and 
total mass (e.g., Colafrancesco 1999; Liang 1999; Feretti 2002, 
2003; Govoni et al. 2001; EnBlin & Rottgering 2002; Cassano 
et al. 2006). These scalings call into question the rarity of halos in 
clusters of low X-ray luminosity, suggesting that the lack of RH 
detections in those clusters may result from the combination of 
the radio power-X-ray luminosity correlation and the sensitivity 
of the radio survey (e.g., Kempner & Sarazin 2001). 

There is also substantial evidence that RHs are found in 
clusters with significant substructure in the X-ray images, 
as well as complex gas temperature distribution, which are 
signatures of cluster mergers (e.g., Schuecker et al. 2001; 
Govoni et al. 2004; Markevitch & Vikhlinin 2001). In particular, 
Buote (2001) provided the first quantitative comparison of the 
dynamical state of clusters with RH, discovering a correlation 
between the RH luminosity at 1.4 GHz and the magnitude of the 
dipole power ratio Pi/ Pq, which is a measure of the cluster’s 
X-ray morphological disturbance. However, these observational 
claims were based on collections of data from the literature and 
not on statistical samples of galaxy clusters. 

An important step was recently obtained through deep ra- 
dio observations of a complete sample of galaxy clusters as 
part of the Giant Metrewave Radio Telescope (GMRT) RH Sur- 
vey (GRHS hereafter; Venturi et al. 2007, 2008). These ob- 
servations confirmed that RH are not ubiquitous in clusters. 
They are found only in ~30% of the X-ray luminous systems 
(Lx(0. 1-2.4 keV) ) 5 x 10 44 erg s -1 ). The sensitivity reached 
by these observations allowed for the first time to place deep 
upper limits on the diffuse radio flux of clusters without giant 
RH and to show that clusters branch into two populations: RHs 
trace the correlation between Pi .4 and Lx, while the upper limits 
on the radio luminosity of clusters with no RH lie about one 
order of magnitude below that correlation (e.g., Brunetti et al. 
2007, 2009). Using several methods to characterize cluster sub- 
structures, it was also shown that clusters with and without RH 
can be quantitatively differentiated in terms of their dynamical 
properties, with RHs always associated with dynamically dis- 
turbed clusters while clusters without RHs are more “relaxed” 
(Cassano et al. 2010). 

Sensitivity is critical in these studies. Indeed, analyses based 
on all-sky surveys, such as the NVSS and WENSS that have a 
sensitivity four to five times worse than the GRHS, do not allow 
to recover a bimodal behavior in the radio-X-ray diagram (e.g., 
Rudnick & Lemmerman 2009 for the WENSS). On the other 
hand, evidence for a bimodal behavior of clusters was recently 
found through a stacking analysis of clusters in the SUMSS 
(Brown et al. 2011). 

More recently, Basu (2012) cross-correlated the Planck ESZ 
cluster catalog (Planck Collaboration et al. 2011a) with radio 
data from the GRHS. He found a correlation between Pi 4 and 
the integrated Sunyaev-Zel’dovich (SZ) effect measurements, 
but did not find a strong indication for a bimodal split between 
RH and radio-quiet clusters. To explain this apparent lack of bi- 
modality in SZ, Basu (2012) suggested that X-ray observations 
could be biased towards the detection of low-mass cool-core 
clusters, whereas SZ selection picks up the most massive sys- 
tems, irrespective of their dynamical states. 

In this paper, we improve on the previous statistical studies 
on the distribution of clusters in the Pi a-L x diagram by using 
RH and clusters with radio upper limits from the GRHS and its 
extension and, when needed, including RHs from the literature. 
Contrary to previous analyses that used cluster X-ray and radio 
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halo luminosities from the literature, we reevaluate the radio 
and X-ray luminosities in a homogeneous way. In particular, 
we derive the X-ray luminosity within /Goo ' 1 and include 
the correction due to the contribution of the cool core (when 
present). Furthermore, since the integrated SZ signal is a more 
robust indicator of the cluster mass than the X-ray luminosity 
(e.g.. Mod et al. 2005; Nagai 2006), we cross-checked our 
sample with the recent Planck SZ (PSZ) cluster catalog (Planck 
Collaboration et al. 2013b) and derived the distribution of 
clusters in the radio-SZ plane. 

In Section 2, we describe the cluster sample. In Section 3, we 
summarize the procedure to derive different cluster quantities 
(radio-halo power, X-ray luminosity, SZ flux, morphological 
parameters), identify cool-core clusters in the sample, and fit 
the scaling relations. In Section 4, we report on the expected 
theoretical scalings. We derive the distribution of clusters in the 
radio-X-ray diagrams in Section 5, and in the radio-SZ (mass) 
diagrams in Section 6. Finally, in Section 7, we give a summary 
and conclusions. 

A ACDM cosmology ( H„ = 70kms _1 Mpc -1 , Q„, = 0.3, 
Q a = 0.7) is adopted. 

2. THE SAMPLE 

The GRHS is a deep, pointed radio survey of clusters se- 
lected from the ROSAT- ESO Flux Limited X-ray (REFLEX; 
Bohringer et al. 2004) and extended ROSAT Brightest Cluster 
Sample (eBCS; Ebeling et al. 1998, 2000) catalogs. These two 
catalogs have almost the same flux limit in the 0. 1-2.4 keV 
band (>3 x 10~ 12 ergs _1 cm -2 ) and their combination yields 
a homogeneous, flux-limited sample of clusters. The GRHS 
consists of 50 galaxy clusters with z = 0.2-0. 4, X-ray lumi- 
nosity Lx > 5 x 10 44 erg s -1 , and declination <5 G 30° for the 
REFLEX sample and 15° ^ S ^ 60° for the eBCS sample. 
With the above selection criteria, the sample is X-ray luminosity- 
limited up to z — 0.25 and X-ray flux-limited at higher redshift 1 2 
(see Figures 1 and 2 in Cassano et al. 2008). 

Recently, we have undertaken an extension of the GRHS by 
considering all clusters in the REFLEX and eBCs catalogs with 
S > —30° with the same z and Lx selection (Kale et al. 2013). 
This extension yields a final sample of 67 galaxy clusters, which 
we refer to as the extended GMRT RH Survey (EGRHS). For 
all clusters in the EGRHS with the radio data already available, 
we searched the ROSAT and Chandra archive and found data 
for a subsample of 40 galaxy clusters: 

1 . 29 with radio upper limits; 

2. 8 with giant radio halos; 

3. 3 with radio halos with ultra-steep spectra. 1 ’ 

In addition to clusters belonging to the EGRHS, we also 
searched in the X-ray archive and found data for 14 clusters 
with RHs from the literature: 

1 . 11 with giant radio halos; 

2. 3 with radio halos with ultra-steep spectra. 

The total sample with radio and X-ray information consists 
of 54 galaxy clusters, whose main properties are reported in 
Table 1: 


1 1 /Goo is the radius corresponding to a total density contrast 500p c (4), where 
p c (z) is the critical density of the universe at the cluster redshift. 

12 This implies a minimum Lx ~ 10 45 erg s' 1 at z ~ 0.35. 

13 We do not include a fourth one, RXCJ15 14.9-1523 (Giacintucci et al. 
2011), which was only recently observed by Chandra , the data analysis is still 
ongoing (S. Giacintucci et al., in preparation). 
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Table 1 

Cluster’s Properties 


Cluster Name 

R-A.J2000 

Decl.j 2000 

z 

£500 

£500, cor 

Lcore/ L 5 OO 

A.4 

X-Ray 

sz 

Upper limits (EGRHS) 

A2697 

00 03 11.8 

-06 05 10 

0.232 

7.29 ± 0.41 

7.29 ± 0.41 

0.34 

<0.41” 

19(H) 

V 

A141 

01 05 34.8 

-24 39 17 

0.230 

6.82 ± 0.27 

6.82 ± 0.27 

0.13 

<0.43" 

33 (H) 

V 

A3088 

03 07 04.1 

-28 40 14 

0.254 

6.97 ± 0.09 

5.63 ± 0.08 

0.46 

<0.43" 

19(C) 

V 

RXCJ0437. 1+0043 

04 37 10.1 

+ 00 43 38 

0.284 

6.99 ± 0.08 

6.15 ± 0.08 

0.45 

<0.65* 

30 (C) 


RXCJ1 115.8 + 0129 

11 15 54.0 

+ 01 29 44 

0.350 

12.69 ±0.11 

8.21 ± 0.10 

0.63“ 

<0.47" 

39 (C) 

V 

A2485 

22 48 32.9 

-16 06 23 

0.247 

3.27 ± 0.07 

3.07 ± 0.07 

0.39 

<0.47* 

20 (C) 


A2631 

23 37 40.6 

+ 00 16 36 

0.278 

8.62 ± 0.70 

8.62 ± 0.70 

0.21 

<0.41” 

15(H) 

V 

A2645 

23 41 16.8 

-09 01 39 

0.251 

4.13 ±0.4 

4.13 ±0.4 

0.43 

<0.59 ? 

35 (H) 

V 

A2667 

23 5140.7 

-26 05 01 

0.226 

12.50 ±0.4 

10.94 ± 0.4 

0.46“ 

<0.45" 

21(H) 

V 

Z348 

01 06 50.3 

+ 01 03 17 

0.255 

6.30 ± 0.60 

4.26 ± 0.60 

0.54“ 

<0.65* 

13(H) 


RXJ0142.0 + 2131 

01 42 03.1 

+ 21 30 39 

0.280 

6.00 ± 0.10 

6.00 ± 0.10 

0.34 

<0.45* 

20 (C) 

V 

A267 

01 52 52.2 

+ 01 02 46 

0.230 

6.29 ± 0.44 

5.94 ± 0.44 

0.36 

<0.34* 

16(H) 

V 

RXJ0439.0 + 0715 

04 39 01.2 

+ 07 15 36 

0.244 

8.05 ± 0.59 

7.69 ± 0.58 

0.36 

<0.46* 

19(H) 

V 

RXJ0439.0 + 0520 

04 39 02.2 

+ 05 20 43 

0.208 

5.35 ± 0.47 

4.05 ± 0.46 

0.55“ 

<0.32* 

12(H) 


A61 1 

08 00 58.1 

+ 36 04 41 

0.288 

4.96 ± 0.64 

4.96 ± 0.64 

0.46 

<0.43" 

17(H) 

V 

Z2089 

09 00 45.9 

+ 20 55 13 

0.235 

5.28 ± 0.41 

3.98 ± 0.41 

0.52“ 

<0.26" 

17(H) 


A781 

09 20 23.2 

+ 30 26 15 

0.298 

5.44 ±0.14 

5.44 ± 0.14 

0.12 

<0.36” 

10(C) 

V 

Z2701 

09 52 55.3 

+ 51 52 52 

0.214 

4.72 ± 0.43 

3.38 ± 0.42 

0.55“ 

<0.35" 

10(H) 


A 1423 

11 57 22.5 

+ 33 39 18 

0.213 

5.35 ± 0.37 

4.76 ± 0.38 

0.41 

<0.38” 

19(H) 

V 

A1576 

12 36 49.1 

+ 63 11 30 

0.30 

6.68 ±0.14 

6.38 ±0.14 

0.24 

<0.64* 

17 (P) 

V 

RXJ 1532.9+3021 

15 32 54.2 

+ 30 21 11 

0.345 

17.94 ± 0.27 

12.28 ± 0.22 

0.64“ 

<0.66” 

22 (C) 


A2146 

15 56 04.7 

+ 66 20 24 

0.234 

5.69 ± 0.04 

5.69 ± 0.04 

0.42 

<0.39'’ 

65 (C) 

V 

A2261 

17 22 28.3 

+ 32 09 13 

0.224 

8.98 ±0.38 

7.79 ± 0.37 

0.43 

<0.32* 

30 (H) 

V 

RXJ2228. 6+2037 

22 28 34.4 

+ 20 36 47 

0.418 

11.71 ±0.20 

11.71 ± 0.20 

0.29 

<0.95" 

20 (C) 

V 

A2537 

23 08 23.3 

-02 11 31 

0.297 

5.48 ±0.14 

4.54 ± 0.07 

0.45 

<0.51" 

39 (C) 

V 

RXJ0027.6 + 2616 

00 27 49.8 

+ 26 16 26 

0.365 

3.52 ±0.11 

3.52 ±0.11 

0.24 

<0.74" 

22 (C) 


Z5699 

13 06 00.0 

+ 26 30 58 

0.306 

4.74 ± 0.08 

4.74 ± 0.08 

0.18 

<0.59" 

26 (C) 


Z5768 

13 11 31.5 

+ 22 00 05 

0.266 

1.66 ± 0.06 

1.66 ± 0.06 

0.10 

<0.36" 

27 (C) 


S780 

14 59 29.3 

-18 11 13 

0.236 

8.68 ±0.10 

6.32 ± 0.05 

0.43 

<0.38” 

40 (C) 

V 

Radio Halos (EGRHS) 

A2744 

00 14 18.8 

-30 23 00 

0.307 

14.73 ± 0.24 

14.73 ± 0.24 

0.17 

18.62 ±0.94 a 

14 (P) 

V 

A0209 

01 31 53.0 

-13 36 34 

0.206 

7.62 ± 0.48 

7.62 ± 0.48 

0.31 

1.99 ±0.21 a 

11(H) 

V 

A2163 

16 15 46.9 

-06 08 45 

0.203 

21.95 ±0.33 

21.95 ±0.33 

0.25 

22.91 ± 1.16 a 

7 (P) 

V 

RXCJ2003.5— 2323 

20 03 30.4 

-23 23 05 

0.317 

9.17 ± 0.09 

9.17 ±0.09 

0.09 

10.71 ± 1.73 b 

50 (C) 

V 

A520 

04 54 19.0 

+ 02 56 49 

0.203 

7.81 ±0.21 

7.81 ± 0.21 

0.18 

2.45 ± 0.1 8 a 

5 (P) 

V 

A773 

09 17 59.4 

+ 51 42 23 

0.217 

7.30 ± 0.57 

7.30 ± 0.57 

0.35 

1.48 ± 0.16 a 

17(H) 

V 

A1758a 

13 32 32.1 

+ 50 30 37 

0.280 

8.80 ± 0.16 

8.80 ±0.16 

0.18 

5.75 ± 0.98 a 

16 (P) 

V 

A2219 

164021.1 

+ 46 41 16 

0.228 

14.78 ±0.19 

14.78 ±0.19 

0.20 

5.63 ± 0.80 a 

16 (P) 

V 

A0521 17 

04 54 09.1 

-10 14 19 

0.248 

8.28 ± 0.07 

8.28 ± 0.07 

0.08 

1.45 ± 0.13" 

39 (C) 

V 

A697" 

08 42 53.3 

+ 36 20 12 

0.282 

13.04 ±0.61 

13.04 ±0.61 

0.33 

1.51 ±0.14‘ 

28 (H) 

V 

A1300 t/ 

11 31 56.3 

-19 55 37 

0.308 

11.47 ±0.37 

11.47 ±0.37 

0.18 

3.8 ± 1.43P 

9 (P) 

V 

Radio Halos (literature) 

clooi6+i6 

00 18 33.3 

+ 16 26 36 

0.541 

15.54 ± 0.28 

15.54 ± 0.28 

0.16 

5.01 ±0.31 a 

43 (P) 

V 

A1914 

14 26 03.0 

+ 37 49 32 

0.171 

11.17 ±0.13 

10.25 ± 0.13 

0.39 

5.62 ± 0.43 a 

9 (P) 

V 

A665 

08 30 45.2 

+ 65 52 55 

0.182 

8.36 ± 0.09 

8.30 ± 0.07 

0.22 

2.51 ±0.21 a 

38 (P) 

V 

A545 

05 32 20.2 

-11 31 54 

0.154 

6.31 ± 0.09 

6.31 ± 0.09 

0.23 

1.41 ±0.22 a 

14 (P) 

V 

Coma 

12 59 48.7 

+ 27 58 50 

0.023 

3.39 ± 0.03* 



0.76 ± 0.06° 


V 

A2256 

17 03 43.5 

+ 78 43 03 

0.058 

4.44 ± 0.02 

4.44 ± 0.02 

0.22 

0.85 ± 0.08 d 

17 (P) 

V 

Bullet 

06 58 29.2 

-55 57 10 

0.296 

22.54 ± 0.52 

22.54 ± 0.52 

0.22 

23.44 ± 1 .5 l e 

5 (P) 

V 

A2255 

17 1231.0 

+ 64 05 33 

0.081 

3.31 ± 0.03 

3.31 ± 0.03 

0.14 

0.81 ±0.17 f 

15 (P) 

V 

A2319 

19 20 45.3 

+ 43 57 43 

0.056 

7.96 ± 0.08 

7.87 ± 0.08 

0.31 

2.45 ±0.198 

3 (P) 

V 

MACS J07 17.5+3745 

07 17 33.8 

+ 37 45 20 

0.548 

24.05 ± 0.22 

24.05 ± 0.22 

0.17 

52.48 ± 20.56 h 

60 (C) 

V 

A 1995 

14 52 50.4 

+ 58 02 48 

0.319 

6.03 ± 0.08 

6.03 ± 0.08 

0.43 

1.66 ± 0.23 a 

48 (C) 

V 

MACSJ1149.5+2223 (/ 

11 49 34.3 

+ 22 23 42 

0.544 

15.50 ±0.29 

15.50 ± 0.29 

0.18 

2.29 ± 0.95 m 

19(C) 

V 

PLCKG 17 1 .9-40. 7 1 * 

03 12 57.4 

+ 08 22 10 

0.270 

11.28 ±0.02** 



4.90 ± 1.35” 


V 

A754" 

09 08 50.1 

-09 38 12 

0.054 

4.75 ± 0.033 

4.75 ± 0.033 

0.23 

0.63 ± 0.07° 

8(P) 

V 


Notes. The first part of the table contains clusters with radio upper limits belonging to the EGRHS (from Venturi et al. 2008; Kale et al. 2013); the second part contains 
clusters with giant RHs belonging to the EGRHS; and the third part contains clusters with giant RHs not belonging to the EGRHS; clusters marked with u are those hosting 
USSRH (or > 1.5). Columns: (1) cluster name; (2) and (3) cluster right ascension and declination, respectively, in J2000 coordinates, taken from the X-ray catalogs (see 
Section 2); (4) cluster redshift; (5) 0. 1-2.4 keV band cluster X-ray luminosity within R 500 ; (6) 0.1-2.4 keV band cluster X-ray luminosity within R 500 corrected for the 
contribution of the cool-core; (7) the ratio between the X-ray luminosity within the core and the total luminosity within R500, cool-core clusters are indicated with cc \ (8) 
^-corrected radio halo power at 1.4 GHz; (9) X-ray exposure in ks, with P = ROSAT PSPC, H = ROSAT HRI, and C = Chandra ACIS-I; (10) the symbol indicates 
the clusters present in the 15.5 month Planck catalog. References for the radio halo powers are v Venturi et al. (2008); k Kale et al. (2013); q Guglielmino 2008, private 
communication; r Russell et al. (2011); a S. Giacintucci et al. (in preparation); b Giacintucci et al. (2009); c Kim et al. (1990); d Clarke & Ensslin (2006 ); e Liang et al. 
(2000); 7 Govoni et al. (2005); 8 Farnsworth et al. (2013); h Bonafede et al. (2009); van Weeren et al. (2009); ' Dallacasa et al. (2009 ); 1 Macario et al. (2010); m Bonafede 
et al. (2012); " Giacintucci et al. (2013); 0 Macario et al. (2011); p Venturi et al. (2013). * the X-ray luminosity of the Coma cluster is taken from O’Hara et al. (2006); ** 
the X-ray luminosity of PLCKG171.9— 40.7 is taken from Planck Collaboration et al. (2011c). 
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1. 8 RH from the EGRHS; 

2. 29 clusters with radio upper limits from the EGRHS; 

3. 6 RH with ultra-steep spectra (a > 1.5; USSRH hereafter); 
not to be compared to the upper limits (which were scaled 
at 1.4 GHz with a = 1.3); 

4. 2 RH from the literature, A 1995 and the Bullet cluster, 
which are in the same redshift and X-ray luminosity range 
of the EGRHS; 

5. 9 RH from the literature, which do not fulfill the EGRHS 
selection criteria (in redshift and X-ray luminosity). 

The inclusion of RH from the literature is useful to have a 
sufficient leverage in radio/X-ray luminosities that may help to 
derive the scaling relations; however, the comparison between 
halos and upper limits will be performed only for clusters of 
EGRHS, which are in the same redshift range. 

3. DATA ANALYSIS 

In this section, we briefly describe the procedures undertaken 
to derive the radio and X-ray luminosities of clusters, to identify 
cool core clusters, and to analyze the cluster dynamical status. 
We also report measurements of the SZ signal found in the 
literature. 

3.1. Radio Power of Halos 

Table 1 reports the radio halo powers and upper limits for the 
clusters in the sample. We refer to the literature information (see 
notes to the table) for details and for the radio images. 

For 12 RH clusters in the sample we re-analyzed archival 
1.4 GHz VLA-C and VLA-D array data (S. Giacintucci et al., in 
preparation). Another six clusters with giant RHs are published 
by our group and the data analysis was carried out following 
procedures similar to those described in this section. For the 
remaining seven RHs, the radio flux densities were taken from 
the literature. For those clusters we analyzed (12+6). the flux 
density of the radio halo was measured from low angular 
resolution images obtained after subtraction of the contribution 
of the individual radio sources embedded in the diffuse emission. 
In particular, we identified the discrete radio galaxies in (or 
projected onto) and around the cluster region using the higher 
resolution images produced with the C-array datasets (when 
available). The discrete radio sources were subtracted from 
the u- v datasets, and the resulting u-v visibilities were then 
used to image the radio halo emission at low resolution. In 
those cases where high-resolution data were not available to 
evaluate and subtract the contribution of individual sources, we 
produced images gapping the innermost region of the u-v plane 
and used only the remaining long baselines (>1-2A), which 
contain information on structures on angular scales smaller than 
the underlying large-scale radio halo. For extended sources, we 
used sets of images with different resolutions and/or u-v ranges 
to determine their total extent and morphology. For each source 
and each cluster, we carefully checked that the total flux density 
subtracted from the u-v data is consistent with the flux density 
measured on the images. 

We measured the total flux density of the radio halos starting 
from the 3er contour level in the final images, then we progres- 
sively increased the extraction region until the integrated flux 
density reached a maximum value, and considered this maxi- 
mum value as the total flux density of the halo. This procedure 
leads to an average increase of the halo flux density by only 
~5% with respect to the value within the 3er isocontour. 


Finally, we calculated the corresponding radio power at the 
cluster redshift and applied a /correction (1 + z) (1_ “- ) , where the 
spectral index is taken from the literature (references in Table 1), 
or was assumed to be a = 1.3 when not available. The errors 
on the diffuse radio flux density,///, account for the uncertainty 
in the calibration of the absolute flux density scale, the error 
due to the noise in the integration area, and the error due to the 
subtraction of the discrete radio sources in the halo region, as 

o-/„ = 7 (<5cai /h) 2 + Inns ^AWO 2 + °7ub’ (!) 

where 5 ca i is typically of the order of 5%-8%, rms is the noise of 
the map, A/cam is the number of independent beams in the halo 
region, and er su b is the error due to the uncertainty in the source 
subtraction. The term ay H does not account for the uncertainty 
due to the missing short spacings in the u-v coverage of the 
interferometric observations, and this may bias the flux densities 
toward lower values. 

Upper limits to the diffuse radio flux of clusters without giant 
RH were reported in Venturi et al. (2008) and Kale et al. (2013) 
and scaled at 1.4 GHz using a spectral index typical of RHs, 
a = 1.3. 

3.2. X-Ray Luminosities 

To derive the cluster X-ray luminosities we use ROSAT data, 
PSPC preferentially and HRI when PSPC data are not available. 
In those cases where ROSAT data are not available, we use 
Chandra data (see Table 1). We derive the X-ray luminosities 
inside /Z 500 centered on the centroid of the X-ray emission for 
all clusters. To estimate //goo for our clusters, we searched 
the literature for information about the X-ray temperature 
and then applied the relations from Arnaud et al. (2005). We 
derived luminosities in the 0. 1-2.4 keV band in three different 
ways: ( 1 ) the entire cluster emission inside //goo, denoted 
as L 500 ; (2) the emission inside the aperture [0.15-l]//goo, 
denoted as Lgoo.nc; (3) Tgoo.cor* the X-ray luminosity inside 
//oo corrected for the excess emission within 0.15 //goo, due 
to the presence of a cooling core. For each cluster, Tsoo.cor 
is computed by performing a fit with a /-model to the cluster 
emission outside 0.15 //goo, fixing r c = 0.15 //goo (assuming that 
r c ~ r c00 i ~ 0.15 //goo, which is ~ 1 00-200 kpc for our sample) 
and evaluating the contribution of the model inside 0.15 //goo- 
When the model fit underestimates the counts in the core, we 
correct the central region by using the fit to the X-ray brightness 
distribution outside 0. 1 5 //goo- We masked the detected point 
sources after a careful inspection of the cluster Chandra images 
(only in four cases we used //OSATPSPC/HRI data). Following 
Russell et al. (2013), bright central point sources were identified 
and masked using the Chandra 5-8 keV images. The flux in the 
masked regions has been replaced by estimates based on 
the cluster best-fit model for the spatial brightness distribution. 
The obtained values of Lgoo and Lgoo, CO r are reported in Table 1 . 

3.3. Identification of Cool-core Clusters 

In this section, we identify cool-core clusters in our sample 
to investigate possible biases that can be induced by cool-core 
clusters on scaling relations and bimodality. 

As first measurement, we consider the X-ray surface bright- 
ness concentration parameter, defined as the ratio between the 
X-ray luminosity within the core region (L c ore , within 0.15 //goo) 
and Lgoo (e.g., Santos et al. 2008; Cassano et al. 2010, with 
slightly different definitions). In the literature, the concentration 
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Figure 1. Left panel: Gore /Goo versus Kq for all clusters in our sample, for those clusters for which we do not find values of Kq in Cavagnolo et al. (2009), we set 
S'o = 10 keV cm~ 2 (at the boundary of the plot). Clusters without giant RHs and clusters with giant RHs are reported as the black and red points, respectively. Central 
panel: t C ool versus Kq for all clusters with /.core/ Goo > 0.4 and Kq < 130 keV cm~ 2 (blue dashed region in the left panel); Right panel: Gore/Goo versus t C ool for all 
clusters with Gore /Goo > 0.4 and Kq < 130 keV cm -2 . 

(A color version of this figure is available in the online journal.) 


parameter has been used for a first identification of cool-core 
clusters in those cases where a spatially resolved spectroscopic 
analysis was not possible (e.g., in the case of high-redshift clus- 
ters; Santos et al. 2008) and to discriminate between merging 
clusters and more relaxed clusters (e.g., Cassano et al. 2010). A 
large value of this parameter indicates a large probability that 
the object has a cool core. 

The derived values of the ratio L corc / Gooo are reported in 
Table 1 . Here we use the concentration parameter in combination 
with the central entropy ( Kq in keV cm -2 ) and the central 
cooling time (f coo i) to identify clusters with a cool-core. Values 
of K 0 < 50keV cm 2 (dashed vertical line in Figure 1) are 
used to identify cool-core clusters (e.g., Cavagnolo et al. 2009; 
Rossetti et al. 2011). 

We inspected the sample of Cavagnolo et al. (2009) to find 
information about the central entropy for our clusters, and in 
Figure 1 (left panel), we report the distribution of clusters in 
the T com /C wo versus Kq diagram; 13 of our clusters are not 
available in the Cavagnolo et al. sample, and in Figure 1 (left 
panel), they are reported with a value of Kq = lOkeV cnG 2 . 
We find the following. 

1. Clusters with giant RHs (open red dots) have L cok /Lsqq < 
0.4 14 and Kq > 90keV cnG 2 . 

2. Five clusters have Kq < 50keV cm" 2 and L cok /L^qq > 

0.5: RXJ1532. 9+3021. RXCJ1 115.8+0129, Z2089, 

RXJ 0439.0 + 0520, Z2701, and one, A2667, has 

C c ore/C500 = 0.46. 

3. Five clusters have 50keV cm " 2 < K q < 130keV cm -2 
and L cole /Lsoo > 0.4: A611, A2261, A3088, A1423, and 
A2537. 

As expected, clusters with giant RHs can be easily identified 
with merging clusters. To better understand whether the 11 
clusters with L core / Goon > 0.4 and Kq < 130 keV cm 2 are 
cool-core or non cool-core clusters, we searched for information 
in the literature about their central cooling time (Cool) (Figure 1, 
central panel). 15 Clusters with Kq < 50keV cnG 2 have Cool < 
2 Gyr, while the others have Cool > 3 Gyr. Figure 1, right 


14 The only exception is A1995 with /.core/ /-too = 0.43 (see also discussion 
in Section 6.2). 

15 We do not find information about t coo \ in the literature for Z2089. 


panel, also shows that clusters with Kq < 50keV cm~ 2 and 
Cool < 2 Gyr all have L c ore/Csoo > 0.5 (with the exception 
of A2667). Therefore, based on the combination of the three 
indicators, we identify cool-core clusters in our sample as those 
with T CO re/ T 500 > 0.5, i.e., clusters that emit more than 50% of 
their L 500 within their cores . 16 Thus, we can conclude that there 
are seven cool-core clusters 17 in our sample (these are marked 
with a cc symbol in Column 7 of Table 1 ). 

3.4. Cluster Sunyaev-Zel’dovich Measurements 

Observations of clusters through their SZ-effect offer a valid 
alternative to X-rays for the measure of the cluster mass, since 
the magnitude of the SZ effect is proportional to the integral 
along the line of sight of the cluster pressure, and hence is 
proportional to the cluster mass. The total SZ signal can be 
defined as 

*a c = D 2 a T SZ = {^~) ( PdV oc 

m e c l J R<t.R &c 

x M Wds T e = f gds M tot T e , (2) 

where D,\ is the angular diameter distance to the system, ay is 
the Thomson cross-section, c is the light speed, m e is the electron 
rest mass, P = n e kT e is the electron pressure, / gas is the gas 
mass fraction, and M tot is the total cluster mass. The integral in 
Equation (2) is performed over a sphere of radius, R\ c , which is 
the radius corresponding to a density contrast, A e p c (z). When 
the integration is performed over a sphere of radius R 500 , the SZ 
signal is denoted with F 500 , which in the following paper will 
have the unit dimension of Mpc 2 . 

For all clusters in Table 1 , we search for information about the 
SZ signal in the recent all-sky PSZ cluster catalog, which con- 
tains all validated clusters from the first 15.5 months of Planck 
satellite observations (Planck Collaboration et al. 2013b). Con- 
sidering only clusters belonging to the EGRHS subsample, we 
find that 11/11 RH clusters and 19 out of 29 clusters with upper 


16 We also consider A2667 as a cool-core cluster, since it has an estimated 
central entropy of Kq m 19 keV cm~ 2 and a central cooling time f coo i « 1 Gyr 
and it is classified as a cool core cluster by Zhang et al. (2007). 

17 The cluster Z348 has no information about Kq in the literature, but since it 
has Gore/Goo = 0.54, we can identify it as a cool-core cluster. 
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limits are contained in the PSZ catalog. Among the 10 clusters 
not present in the PSZ catalog, 5 are cool-core clusters; there- 
fore, only 2 out of 7 cool-core clusters of our sample are detected 
by Planck. The remaining 14 RH clusters from the literature are 
also contained in PSZ catalog. 

We obtain a subsample of 44 clusters (25 halos and 19 upper 
limits) for which Planck measurements of i^oo are available (see 
Table 2). 

For the same clusters, we also find information in the PSZ 
catalog about the values of M5 00 ■ These are obtained from T500 as 
described in Planck Collaboration et al. (2013b, Section 7.2.2) 
and are reported in Table 2. 

3.5. Cluster Dynamical Status 

For clusters belonging to the EGRHS with information about 
T 500 and M 500 (see Table 2) we make use of Chandra archival 
data to determine the cluster dynamical status. 18 We produce 
X-ray images in a standard manner using CIAO 4.3 (with 
calibration hies from CALDB 4.4.1) in the 0.5-2 keV band. 
We adopted an algorithm for an automatic detection of the point 
sources, which are then removed from the images. Following 
Cassano et al. (2010), we study the cluster substructure on the 
RH scale analyzing the surface brightness inside an aperture 
radius of 500 kpc, since we are interested in the cluster 
dynamical properties on the scales where the energy is most 
likely dissipated. We use two methods: the emission centroid 
shift (e.g., Mohr et al. 1993; Poole et al. 2006; O’Hara et al. 
2006; Ventimiglia et al. 2008; Maughan et al. 2008; Bohringer 
et al. 2010) and the surface brightness concentration parameter 
(e.g., Santos et al. 2008). 

The centroid shift, w, is computed in a series of circular 
apertures centered on the cluster X-ray peak and is defined as 
the standard deviation of the projected separation between the 
peak and the centroid in units of R dp as (Poole et al. 2006; 
Maughan et al. 2008) 

r 1 v—' n 1 / 2 1 

” = [^tE<A,-(A)) j ] X-, (3) 

where A, is the distance between the X-ray peak and the centroid 
of the zth aperture. 

Following Santos et al. (2008), the concentration parameter, 
c, is defined as the ratio of the peak over the ambient surface 
brightness, S, as 

S(r < 100 kpc) 

c = . (4) 

5(<500kpc) 

We use the concentration parameter to differentiate galaxy 
clusters with a compact core (i.e., core not disrupted from a 
recent merger event) from clusters with a spread distribution 
of gas in the core (i.e., core disturbed from a recent merger 
episode). 

Cassano et al. (2010) showed that, considering the median 
value of each parameter, w = 0.012 and c = 0.2, it was 
possible to separate the sample between RH merging clusters 
( w > 0.012 and c < 0.2) and more relaxed clusters without 
RHs (w < 0.012 and c > 0.2). We will use these values as 
reference for our sample. 


18 With the exception of A2697 for which Chandra data are not available. 


Table 2 

Observed Cluster SZ Properties 


Cluster Name 

Index 

iog(y 5 oo) 

log(M 50 o) 



(Mpc 2 ) 

(Mo) 


Upper limits (EGRHS) 








A2697 

315 

-4.150 

± 

0.077 

14.78 

± 

0.04 

A141 

599 

-4.379 

± 

0.120 

14.65 

± 

0.07 

A3088 

744 

-4.062 

± 

0.065 

14.83 

± 

0.04 

RXCJ1 115.8 + 0129 

881 

-4.087 

± 

0.087 

14.80 

± 

0.05 

A2631 

297 

-4.029 

± 

0.067 

14.84 

± 

0.04 

A2645 

254 

-4.288 

± 

0.099 

14.70 

± 

0.06 

A2667 

94 

-4.054 

± 

0.055 

14.83 

± 

0.03 

RXJ0142. 0 + 2131 

500 

-4.134 

± 

0.102 

14.78 

± 

0.06 

A267 

541 

-4.301 

± 

0.108 

14.69 

± 

0.06 

RXJ0439.0 + 0715 

640 

-4.181 

± 

0.096 

14.76 

± 

0.05 

A611 

623 

-4.162 

± 

0.081 

14.77 

± 

0.05 

A781 

654 

-4.097 

± 

0.072 

14.80 

± 

0.04 

A1423 

610 

-4.143 

± 

0.064 

14.78 

± 

0.04 

A1576 

460 

-4.143 

± 

0.063 

14.78 

± 

0.04 

A2146 

359 

-4.495 

± 

0.080 

14.58 

± 

0.05 

A2261 

174 

-3.991 

± 

0.048 

14.87 

± 

0.03 

RXJ2228.6 + 2037 

275 

-3.917 

± 

0.072 

14.89 

± 

0.04 

A2537 

247 

-4.120 

± 

0.080 

14.79 

± 

0.04 

S780 

1185 

-3.957 

± 

0.062 

14.89 

± 

0.03 

S780 

1185 

-3.957 

± 

0.062 

14.89 

± 

0.03 

Radio Halos (EGRHS) 








A2744 

26 

-3.778 

± 

0.041 

14.98 

± 

0.02 

A0209 

558 

-3.916 

± 

0.041 

14.91 

dfc 

0.02 

A2163 

19 

-3.374 

± 

0.019 

15.22 

± 

0.01 

RXCJ2003.5— 2323 

46 

-3.967 

± 

0.068 

14.87 

± 

0.04 

A520 

655 

-4.030 

± 

0.062 

14.85 

± 

0.04 

A773 

578 

-4.026 

± 

0.049 

14.85 

± 

0.03 

A1758a 

389 

-3.922 

± 

0.044 

14.90 

± 

0.03 

A2219 

242 

-3.681 

± 

0.026 

15.04 

± 

0.01 

A521 17 

688 

-4.040 

± 

0.070 

14.83 

± 

0.04 

A691 u 

628 

-3.640 

± 

0.032 

15.06 

± 

0.02 

AU00 u 

960 

-3.839 

± 

0.053 

14.95 

± 

0.03 

Radio Halos (literature) 








CL0016+ 1609 

408 

-3.813 

± 

0.077 

14.94 

± 

0.04 

A1914 

224 

-4.045 

± 

0.039 

14.84 

± 

0.02 

A665 

533 

-3.914 

± 

0.037 

14.92 

± 

0.02 

A545 

707 

-4.397 

± 

0.112 

14.64 

± 

0.06 

Coma 

187 

-4.281 

± 

0.030 

14.72 

± 

0.02 

A2256 

407 

-4.135 

± 

0.022 

14.80 

± 

0.01 

A2255 

325 

-4.288 

± 

0.028 

14.71 

± 

0.02 

A2319 

252 

-3.900 

± 

0.020 

14.93 

± 

0.01 

MCSJ0717.5 + 3745 

608 

-3.612 

± 

0.049 

15.05 

± 

0.03 

Bullet 

920 

-3.577 

± 

0.025 

15.09 

± 

0.02 

A 1995 

337 

-4.257 

± 

0.075 

14.71 

± 

0.04 

MCSJ1 149.5 + 2223 u 

765 

-3.824 

± 

0.072 

14.93 

± 

0.04 

PLCKG171.9-40.7 17 

591 

-3.666 

± 

0.039 

15.05 

± 

0.02 

A754 u 

801 

-4.095 

± 

0.023 

14.82 

± 

0.01 

Notes. Columns: (1) radio properties; (2) cluster name; (3) index indicating 
the position in the Planck validation catalog; (4) logarithmic value of F 500 in 
Mpc 2 , with 68% errors; (5) logarithmic value of M 500 in solar masses, with 
68% errors. Clusters marked with u are those hosting USSRHs. All the M 500 
and I 500 values refer to Planck Collaboration et al. (2013b; from the Web 
site: http://szcluster-db.ias.u-psud.fr); with the exception of the cluster PLCK 
G171.9— 40.7, whose values are taken from Planck Collaboration et al. (2011c). 


3.6. Fitting Procedure 

Here we describe the procedure used in the next sections to 
investigate the presence of scaling relations between indepen- 
dent measurements, i.e., the RH power and the cluster ther- 
mal quantities (L 500 , Tsoo.cor, F 500 , and M 5 oo). For each set of 
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observables, we fit a power-law relation using linear regression 
in the log-log space by adopting the BCES-bisector and the 
BCES-orthogonal regression algorithms (Akritas & Bershady 
1996) which treat the variables symmetrically and take into ac- 
count measurement errors in both variables and intrinsic scatter 
in the data. Performing Monte Carlo simulations to test the per- 
formances of different regression methods, Isobe et al. (1990) 
recommended the use of BCES-bisector in the case one would 
like to treat the variables symmetrically. Consequently, we will 
consider the BCES-bisector as the reference method. 

Since we also have upper limits on P\a, in those cases where 
upper limits and detections are not clearly separated, we also use 
a regression analysis based on the parametric EM (Expectation- 
Maximization) algorithm that is implemented in the ASURV 
package (Isobe et al. 1986) and deals with “censored data,” 
upper limits. 

Assuming a linear relation of the form Y = aX + b, and a 
sample of N data points (T;, Xj) with errors a Yl and o Xl , we 
estimate the raw scatter using the error weighted orthogonal 
distances to the regression line (e.g., Pratt et al. 2009; Biffi et al. 
2013): 

1 N 

°L = E u, ‘ (y '- - aX ‘ ~ b ?’ (5) 

where 

= 1M 2 

' (W)£f = i iM 2 

Since we are dealing with a limited sample, the regression 
line obtained for our data is a sample regression line that can 
deviate from the (unknown) true regression line. To evaluate the 
variation of our best-fit relation about the true regression line, we 
estimate the 95% confidence interval for the mean value of ( Y ) 
at a given X, i.e., the area that has a 95% chance of containing 
the true regression line. For a given value of the X variable the 
95% confidence region around the mean (Y) (which is given by 
the best-fit relation: (Y) = aX + b) is (Y) ± AT, where 


and er 2 = Oy. + a 2 o\. . (6) 


yy (Y, ~ Y m ) 2 / 1 (X - X m ) 2 \ 

U N ~ 2 '\\ N + T,ti( x i-x m ) 2 /’ 

(7) 

where for each observed X,, Y m = aX, + b, and X m = 

Lilt X i/N- 

4. EXPECTED SCALING RELATIONS 

Scaling relations between the synchrotron radio power of ha- 
los and the cluster thermal properties (mass. X-ray luminosity, 
temperature) are expected in theoretical models for the forma- 
tion of giant RHs. In this section, we briefly summarize the basic 
theoretical expectations for the scalings. 

4.1. Secondary Models 

In the simplest scenario for the formation of giant RHs in 
clusters, the electrons responsible for the synchrotron emission 
are secondary products of the hadronic interaction between 
thermal and cosmic ray protons. In this model, following the 
formalism by Kushnir et al. (2009), the scaling between the 
synchrotron radio power and the cluster [0. 1—2.4] keV X-ray 

a^+0.5 

luminosity is expected to be v P„ syn oc L“y , where a i is 


AT = ±1.96 

\ 


the slope of the L-T relation. For a L — 2.5-3 (e.g., Markevitch 
1998; Arnaud & Evrard 1999; Reiprich & Bohringer 2002; Pratt 
et al. 2009), one obtains 

v P„ syn oc Z^ 58 " 1 ' 46 . (8) 

This is valid under the assumption that the relevant radiation 
losses for the secondary electrons are synchrotron losses, i.e., 
assuming that the average magnetic field strength in the halo 
volume is B > Bcmb — 3.2(1 + z) 2 /iG. Lower magnetic field 
values are disfavored by the combination of Planck and Fermi 
data with radio observations (e.g., Jeltema & Profumo 2011; 
Planck Collaboration et al. 2013a; Brunetti et al. 2012). 

Since T500 is found to scale as T500 oc £^ 02±007 (e.g., Planck 
Collaboration et al. 2011b), this model predicts 

v P„ syn oc (T500) 1 ' 55-1 ' 43 . (9) 

4.2. Turbulent Re-acceleration Models 

In the case of the turbulent re-acceleration scenario, the 
derivation of scaling relations is less straightforward, due to 
our poor knowledge of the details of the microphysics of the 
ICM. A simple approach to derive scaling relations in this 
model is presented in Cassano et al. (2007). Under quasi- 
stationary conditions, the energy flux of the turbulence which 
goes into relativistic electrons is reradiated via synchrotron and 
IC mechanisms. The injection rate of the turbulence generated 
during a merger in the RH volume can be estimated as s, oc 
pH x v 2 /r cms , where pH is the ICM mean density in the RH 
volume, v is the cluster-cluster impact velocity, and x cros is 
the cluster crossing time. As in the case of secondary models, 
it is assumed that the ratio between the energy densities in 
relativistic particles and thermal plasma does not change in any 
systematic way with cluster mass (or temperature) among RH 
clusters. Under this hypothesis, the synchrotron radio power 
is v Py yn oc (M H (Jh)/IF(z, M h , b H ), where ,F(z, M H , b H ) = 
[1 +(3.2(1 + z) 2 / B h) 2 ], and Mh, cth, and Bh are the total cluster 
mass, the cluster velocity dispersion, and the average magnetic 
field strength within the RH size (Rh), respectively (Cassano 
et al. 2007). The expression T is constant in the asymptotic 
limit Bjj 3>> P 2 mb , or when the magnetic field in the RH region 
is independent of the cluster mass. In this case, v P„ syn oc M\f . 
Assuming the scalings Mh oc R 2 H n (Cassano et al. 2007) and 
Rh oc P b Q 0 (Basu 2012), one has 

v P? a cx M500, GO) 

and considering the scaling M500 oc T^ 1 ' 74 , one has 

v P v syn cx Tfo 3 , (11) 

which is steeper than that predicted by “secondary models.” Re- 
acceleration modes also allow the case B 2 H Xg. B 2 mb , without 
tension with y -ray upper limits (e.g., Brunetti et al. 2012), and 
in this case one has T~ x oc M 2 ^” , which implies a correlation 
even steeper than that obtained in the previous case. 

Besides the details of the slopes of the thermal-non-thermal 
scaling relations expected from a different origin of the emitting 
electrons, an important difference between the two scenarios 
is the expected dispersion of the correlations. Re-acceleration 
models predict a variety of spectral shapes of RHs, including 
very steep spectra (e.g., Cassano et al. 2006; Brunetti et al. 
2008), which imply a substantial dispersion in the correlations 
(Kushnir et al. 2009; Brunetti et al. 2009) and an increase of the 
scatter at low observing frequency (Cassano 2010). 
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Figure 2. Left panel: distribution of clusters in the P\ 4 — L 500 plane. Right panel: distribution of clusters in Pi .4 — Lsoo.cor plane. In both panels, different symbols 
indicate halos belonging to the EGRHS (blue filled dots); halos from the literature (black open dots); halos with very steep spectra (USSRH, green asterisks); A 1995 
and Bullet cluster (blue stars); cool core clusters belonging to the EGRHS (magenta arrows). Best-fit relations to giant RHs only (black lines) and to all RHs (including 
USSRH, green dashed lines) are reported. The 95% confidence regions of the best-fit relations obtained for giant RHs only are also reported (shadowed regions). 

(A color version of this figure is available in the online journal.) 


5. RADIO-X-RAY LUMINOSITY CORRELATION 
AND THE BIMODALITY 

It is well known that the radio luminosity of halos at 1 .4 GHz 
scales with the X-ray luminosity of the hosting clusters (e.g., 
Liang et al. 2000; Feretti 2002, 2003; EnBlin & Rottgering 
2002; Cassano et al. 2006; Brunetti et al. 2009; Giovannini 
et al. 2009). This correlation has been used to claim that a 
correlation should also exist between the radio power and the 
virial mass of the host cluster (e.g., Cassano et al. 2006). Deep 
upper limits to the radio flux density of clusters with no RH 
emission at 610 MHz, which were a factor of ~3-20 below the 
correlation, were obtained from the GRHS and its extension 14 
allowing to validate the correlation itself and to discover the 
radio bimodality (e.g., Brunetti et al. 2007). 

In previous papers, the distribution of galaxy clusters in the 
radio-X-ray luminosity diagram, and the scaling relation be- 
tween the two quantities, were based on non-homogeneous radio 
and X-ray measurements. In particular, the radio luminosities of 
halos were collected from the literature and X-ray luminosities 
were taken from RASS-based cluster catalogs. Here we recom- 
puted the radio flux densities of well known RHs by reanalyzing 
observations from the archives (as outlined in Section 3.1). For 
all clusters, we computed the 0. 1-2.4 keV X-ray luminosities 
within /Coo from pointed ROSAT and Chandra observations (see 
Section 3.2). 

In Figure 2, (left panel) we show the distribution of clusters 
in the Pi 4 — L 500 diagram. We report with different colors 
clusters belonging to the EGRHS (blue points and blue and 
magenta arrows) and halos from the literature (black points). 
This is necessary, since the comparison between RH powers 
and upper limits makes sense only for those clusters observed 


19 Previous attempts to compare upper limits and the correlation can be found 
in Dolag (2006). 


within the same redshift range, and this is possible only for 
clusters belonging to the EGRHS. Halos from the literature 
follow the same distribution of halos from the EGRHS, and thus 
we use them to draw the correlation. RH clusters appear to follow 
a well-defined correlation between the halo radio power and 
Z. 500 - Being steeper than other halos, ultra-steep spectrum RH 
(green asterisks) are, in general, under-luminous with respect 
to this correlation. Remember that the position of USSRH in 
the P l 4 — L 500 diagram cannot be compared with that of the 
upper limits as the latter were scaled at 1.4 GHz using a — 1.3. 
We find a bimodal distribution of clusters with the presence of 
two distinct populations, that of radio-halo clusters and that of 
radio-quiet clusters. For values of P 500 > 5 x 10 44 erg s _1 , 
clusters with upper limits to the radio power (blue and magenta 
arrows) are all located below the 95% confidence region of the 
correlation. 

As the EGRHS is based on X-ray-selected clusters, one may 
suspect that the bimodality could be caused by the presence 
of cool-core clusters, which are brighter in X-ray and do not 
host giant radio-halos. With the idea to test the bimodality 
against the presence of cool-core clusters in the EGRHS, 
we derive the distribution of clusters in the Pi 4 — Lsoo.cor 
diagram (Figure 2, right panel). We highlight the position of 
cool-core clusters (identified as outlined in Section 3.3, magenta 
arrows in Figures 2). As expected, the X-ray luminosity of 
cool-core clusters is significantly reduced going from Z .500 

tO p500,cor- 

However, the bimodal behavior in the halo radio power also 
remains in the Pi .4 — Psoo.cor diagram. Also, in this case, if we 
restrict to clusters with P 5 oo.cor ^ 5 x 10 44 erg s _1 , upper limits 
are all below the 95% confidence region of the correlation. We 
may thus conclude that the observed radio bimodality is not 
driven by the presence of cool-core clusters without diffuse 
radio emission in the EGRHS. We fit the observed Pi .4 — P 500 
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and Pi 4 — L goo, cor relation with a power-law of the generic form 

log ( „ . Pl A T ) = B log ( ,. Lx t)+A, (12) 

\ io 24 - 5 Watt Hz~' / “ V 10 45 erg s _1 / 

where Lx is L 500 or Psoo.cor- The fit was performed using 
linear regression in the log-log space by adopting both the 
BCES-bisector and BCES-orthogonal methods (as discussed in 
Section 3.6). The results of the fit, together with that from 1000 
bootstrap resamples, are reported in Table 3. The slope of the 
correlation is ~2. 1 ± 0.2 and ~2.2 ± 0.2 in the BCES-bisector 
and BCES-orthogonal cases, respectively, consistent with that 
found in previous studies (e.g., Brunetti et al. 2009). The best- 
fit relation has a lower normalization and a larger a raw when 
USSRH are included in the fit (see Table 3). 

6 . RADIO-SZ SCALING RELATIONS 

As discussed in Section 3.4, observations of clusters through 
their SZ-effect may provide a powerful method to measure 
the cluster masses. Recently, Basu (2012) found a correlation 
between the radio power of clusters with RHs and the integrated 
Compton parameter derived from the Planck ESZ catalog 
(Planck Collaboration et al. 2011a) in the form Pi .4 oc Y$ R , 
where T 5 r 500 is the integral of the SZ signal within a radius of 
5 /Goo . 2 ' 1 Basu (2012) found indication for a weaker or lack of 
bimodality based on the fact that only four clusters from the 
GRHS with radio upper limits were found in the Planck ESZ 
catalog, while almost all RH of the GRHS have counterparts in 
the same catalog. Basu (2012) suggested that a possible reason 
for the lack of bimodality in SZ could be due to the fact that 
X-ray selected cluster samples are biased towards the detection 
of X-ray luminous, but not necessarily massive, clusters, while 
the SZ tends to be more “mass-limited.” In this picture, clusters 
with radio upper limits that are not detected by Planck should 
be less massive systems (with respect to those hosting giant 
RHs) with cool-cores. These clusters would appear brighter 
in X-ray because of the presence of a cool core, causing an 
apparent bimodality in the P\ 4 — L 500 plane. However, as we 
have shown in Section 5, even when we consider the X-ray 
luminosity excising the cool core, we find a clear bimodality in 
the radio-X-ray plane (Figure 2, right panel). 

The all-sky PSZ catalog that we are using in this paper, is six 
times the size of the Planck ESZ catalog (Planck Collaboration 
et al. 2013b) used by Basu (2012), and ~80% complete for 
M 500 > 6 x 10 14 M 0 at z — 0.2-0.35, typical mass and 
redshift ranges of the EGRHS clusters. In Figure 3, we show 
the distribution of the 44 clusters of our sample belonging to 
the PSZ catalog (see Section 3.4), in the Pi 4 — M 500 (left panel) 
and Pi 4 — T 500 (right panel) diagrams. We show with different 
colors clusters belonging to the EGRHS (blue points and blue 
and magenta arrows), halos from the literature (black points), 
and halos with ultra-steep radio spectra (green asterisks). The 
comparison between RHs and upper limits can be performed 
only for clusters belonging to the EGRHS, while the RHs from 
the literature are added to better determine the correlations. 
We find clear correlations between Pi 4 and Mgoo and Tgoo 
parameters. Using the BCES regression method, we fit the 
observed Pi 4 — T 500 and Pi .4 — M 500 relation with the following 
power laws: 

log ( ,,, PlA r )=Plog ( — ^ 500 ) + A (13) 

0 V 10 24 5 Watt Hz- 1 / S V 10 _ 4 Mpc 2 / 

20 F 5 /? 500 can be rescaled to Y 500 for the fiducial GNFW model as 
Y 5 r 500 = 1.79 x Y 500 (Amaud et al. 2010). 


Table 3 

Best-fit Parameters of Scaling Relations 


Method 

B 

err(B) 

A 

err(A) 

°raw 

r s 

P 

Pl A ~ £500 

RH+USS 
BCES bisector 

2.11 

0.20 

0.088 

0.056 

0.23 

0.83 

2.32 x 10- 7 

Bootstrap 

2.11 

0.21 

0.083 

0.058 




BCES orthogonal 

2.35 

0.25 

0.094 

0.058 




Bootstrap 

2.37 

0.31 

0.089 

0.062 




RH only 
BCES bisector 

2.10 

0.17 

0.181 

0.048 

0.20 

0.95 

1.03 x 10- 7 

Bootstrap 

2.11 

0.19 

0.176 

0.049 




BCES orthogonal 

2.20 

0.18 

0.185 

0.049 




Bootstrap 

2.21 

0.23 

0.180 

0.049 






Pi . 4 — P 500 , 

cor 




RH+USS 
BCES bisector 

2.11 

0.20 

0.091 

0.056 

0.23 

0.83 

2.32 x 10~ 7 

Bootstrap 

2.12 

0.22 

0.085 

0.060 




BCES orthogonal 

2.35 

0.25 

0.098 

0.058 




bootstrap 

2.38 

0.31 

0.094 

0.065 




RH only 
BCES-bisector 

2.11 

0.16 

0.186 

0.048 

0.20 

0.95 

1.03 x 10~ 9 

Bootstrap 

2.11 

0.18 

0.184 

0.050 




BCES orthogonal 

2.20 

0.18 

0.190 

0.049 




Bootstrap 

2.22 

0.22 

0.187 

0.052 




Pi .4 — ALoo 

RH+USS 
BCES bisector 

3.70 

0.56 

0.009 

0.074 

0.37 

0.73 

3.98 x 10~ 5 

Bootstrap 

3.73 

0.64 

0.011 

0.079 




BCES orthogonal 

5.05 

0.99 

0.002 

0.094 




Bootstrap 

5.27 

1.33 

-0.002 

0.107 




RH only 
BCES-bisector 

3.77 

0.57 

0.125 

0.076 

0.35 

0.81 

2.50 x 10~ 5 

Bootstrap 

3.84 

0.66 

0.126 

0.079 




BCES orthogonal 

4.51 

0.78 

0.129 

0.087 




Bootstrap 

4.62 

0.90 

0.131 

0.092 




Pi .4 - Uoo 

RH+USS 
BCES bisector 

2.02 

0.28 

-0.131 

0.070 

0.35 

0.74 

2.66 x 10~ 5 

Bootstrap 

2.03 

0.30 

-0.133 

0.069 




BCES orthogonal 

2.48 

0.43 

-0.167 

0.089 




Bootstrap 

2.55 

0.51 

-0.177 

0.100 




RH only 
BCES-bisector 

2.05 

0.28 

-0.014 

0.068 

0.32 

0.83 

1.26 x 10~ 5 

Bootstrap 

2.07 

0.30 

-0.016 

0.072 




BCES orthogonal 

2.28 

0.35 

-0.027 

0.073 




Bootstrap 

2.30 

0.38 

-0.030 

0.079 




RH+UL 
EM algorithm 

2.77 

0.54 

-0.55 

0.13 




RH only 
EM algorithm 

1.70 

0.26 

0.006 

0.068 





Note. The last two columns gives the Spearman’s rank correlation coefficient, 
r s , and the related probability of no correlation. 


and 


log( 


1.4 




10 24 - 5 Watt Hz 


+ A. (14) 


Results of the fits, together with those from 1000 bootstrap 
resamples, are reported in Table 3. The slope of the Pi 4 — Loo 
correlation is close to ~2, consistent with that found by Basu 
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Figure 3. Distribution of clusters in the Pi .4 — M500 (left panel) and in the P1.4 — J500 diagrams (right panel). In both panels, the different symbols are as in Figure 2 . 
Best-fit relations to giant RHs only (black lines) and to all RHs (including USSRH, green dashed lines) are reported. The dashed line in the right panel marks the value 
I500 = 6 x 10 -5 Mpc 2 . 

(A color version of this figure is available in the online journal.) 


(2012); it is 2.05 ± 0.28 when the BCES-bisector method is 
used, and 2.28 ± 0.35 when the BCES-orthogonal method is 
adopted. The slope of the Pi 4 — M 500 correlation is 3.77 ± 0.57 
and 4.51 ± 0.78 in the case of the BCES-bisector and BCES- 
orthogonal methods, respectively, both steeper than previous 
estimates because of the different definitions of the cluster 
masses (within a fixed size of 3 Mpc, Feretti 2003; or the virial 
mass, Cassano et al. 2006). 

At variance with Basu (2012), we find a clear bimodal 
behavior of clusters in both diagrams. For Mgoo > 5.5xl0 14 M G 
and for F 500 > 6 x 10 ~ 5 Mpc 2 , all clusters with radio upper 
limits are well below the 95% confidence region of the best-fit 
correlations. For the sake of completeness, for the Pi 4 — Y$z 
relation, we also performed a regression analysis by making 
use of the parametric EM algorithm that also deals with upper 
limits (see Section 3.6). This allows us to evaluate the effect of 
the radio upper limits on the best-fit correlation, and thus to test 
the reliability of the correlation and the presence of a bimodal 
behavior in the cluster radio powers. The best-fit values are 
reported in Table 3 and the best-fit correlations obtained for 
giant RHs only and for giant RHs plus upper limits are shown in 
Figure 4 (solid and dashed line, respectively) together with the 
95% confidence region of the RH-only correlation. All upper 
limits (with just one exception) lie below the 95% confidence 
region, and the two best-fit relations obtained by considering 
RHs plus upper limits or only RHs differ both in slope and in 
normalization. 

Our statistical analysis suggests two distinct populations of 
clusters: those with giant RHs, occupying the region of the 
correlation, and those without giant RHs, separated from that 
region. 

6.1. N on-detected Planck Clusters in the P \ .4 — F 500 

With the aim of evaluating the possible position, in the 
Pi .4 — F 500 diagram, of EGRHS clusters not contained in the 
15.5 month PSZ validation catalog, we make use of the corre- 
lation between F 5 qo and the core-excluded X-ray luminosity, 



Y 500 [Mpc 2 ] 


Figure 4. Distribution of clusters in P\ .4 — F500 plane. The symbols are as in 
Figure 2 . Best-fit relations to giant RHs (black solid line) and to giant RHs plus 
upper limits (dashed line) are also shown. The shadowed region show the 95 % 
confidence region of the best-fit correlation for giant RHs. 

(A color version of this figure is available in the online journal.) 


T500, nc (Section 3.2). By using the Planck-XMM-Newton 
archive sample, which comprises 62 clusters with the high- 
est quality X-ray and SZ data set currently available (Planck 
Collaboration et al. 2011b), we derive the 0. 1-2.4 keV X-ray 
luminosity between [0. 1 5— 1 ] /?5oo ( 7 - 500 , nc, hereafter) and ob- 
tained the following I 500 — 7 - 5 oo, n c correlation: 


Kz)- 2 / 3 y 5 00 = 


A 


/ h(z) 7/3 7- S00 , nc \ 
v7 x 10 44 ergs -1 / 


Mpc 2 , 


(15) 
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Y 50 o [Mpc 2 ] 


Figure 5 . Comparison between the observed values of F500 (in abscissa) and 
those predicted by the F 500 — £ 500 , nc relation, F 5 oo,pred (in ordinate) for RH 
clusters (red points), and clusters with radio upper limits (black points). The 
black solid line shows the one-to-one trend. 

(A color version of this figure is available in the online journal.) 

where h(z) = VO ffl (l + z) 3 + O a , A = i()- 3 - 795±0014 , and 
B — 1.094 ± 0.039. We thus derive L 500 .no for all clusters 
in Table l 21 and then apply Equation (15) to estimate their 
T 500 parameters. To test the consistency of this approach, we 
compare the “observed” and “predicted” values of T 500 for the 
clusters present in the PSZ catalog. Such a comparison is shown 
in Figure 5; the data are consistent with a one-to-one trend (with 
increasing scatter at lower values), suggesting that indeed we 
can apply this procedure to get reliable estimates of the T 500 for 
clusters not contained in the PSZ catalog. 

In Figure 6 , we show the distribution of all clusters in 
the L\ 4 - Z 500 diagram, including those that are actually not 
observed by Planck (dashed arrows). As expected, the bulk of 
clusters missing in the PSZ catalog is in the region of lower T 500 
values and with M 500 < 5.5 x 10 14 M Q , where the PSZ catalog 
is only marginally complete (the completeness is ~20%). There 
are, however, two exceptions : RXCJ 1532.9+3021, a luminous 
cool core cluster, and RXCJ043 7. 1+0043, which are expected 
in the region of massive clusters. 

6.2. On the P\a — T 500 Scaling Relation 

If we focus on clusters belonging to the EGRHS and also 
consider the two clusters (A 1995 and the Bullet cluster) which 
are in the same X-ray luminosity and redshift range of the 
EGRHS clusters , 22 we find a segregation of clusters in the 
P\ 4 — T 50 o diagram for T 500 > 6 x 1 0 2 Mpc 2 : clusters 
with RHs follow a trend between their radio power and the 
cluster SZ parameter, while clusters without RHs populate the 
region of radio upper limits, which is a factor of ~5-7 below 


21 With the exceptions of Coma, MACS 1149.5+2223, and PLK1 7 1.9-40.7, for 
which the L 500 were taken from the literature (see Table 1). 

22 A 1995 belongs to the NORAS survey which has a slightly lower flux limit 
with respect to the eBCS used to select the GRHS; the Bullet cluster is in the 
south and not easily accessible for the GMRT. 



Y soo [Mpc 2 ] 


Figure 6 . Distribution of clusters in the plane P\a — Fsoo- The symbols are as in 
Figure 2, with the dashed arrows indicating the predicted positions of clusters 
currently not present in the Planck catalog. The best-fit to giant RHs only (black 
solid line) and to giant RHs plus USSRH (green line) are also shown. The 
shadowed region show the 95% confidence region of the best-fit correlation for 
giant RHs only. 

(A color version of this figure is available in the online journal.) 


the correlation (Figure 3, right panel). On the other hand, for 
T 500 < 6 x 10 ~ 5 Mpc 2 , upper limits are not deep enough and lie 
within the 95% confident region of the best-fit correlation. 

In order to better understand this behavior of clusters and 
shed light on the mechanism responsible for the formation of 
giant RHs in clusters, we looked at the dynamical properties of 
clusters in the P\ 4 — T 500 diagram, adopting the centroid shift 
variance, w, and the surface brightness concentration parameter, 
c, to differentiate between merging (w < 0.012 and c > 0.2) 
and more relaxed ( w > 0.012 and c < 0.2) systems (see 
Section 3.5). 

For T 500 < 6 x 10 ~ 5 Mpc 2 , we find that four clusters with 
radio upper limits detected by Planck are merging clusters 
(Figure 7, left panel). Unfortunately, the L 1.4 — I 500 correlation 
predicts an RH power for those clusters that is below the 
sensitivity of the current radio data. In light of our results, 
the apparent lack of a giant RH in the merging but relatively 
low-mass cluster A2146 (Russell et al. 2011) is not surprising, 
because even if a halo is present in this cluster, it may not 
be luminous enough to be detected. The only RH cluster with 
T 500 < 6 x 10 2 Mpc 2 is A1995, which in Figure 7 (left 
panel) is located in a region generally populated by “relaxed” 
clusters. However, A1995 is a merging system, but the merger is 
happening mainly along the line of sight (Boschin et al. 2012), 
and for this reason, its position in the c—w diagram is likely 
biased by projection effects. 

Clusters with T 500 > 6 x 1 0“ s Mpc 2 show a clear segregation 
in their dynamical properties. All clusters with detected giant 
RHs are clearly merging systems, while the majority of clusters 
with upper limits (~80%) are more relaxed (Figure 7, right 
panel). The presence of a segregation in the dynamical state 
of clusters with detected and non-detected RH strengthens the 
separation of clusters in the L\ 4 — T 500 diagram and suggests 
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Figure 7. Distribution of clusters in the plane c—w. Clusters with Jsoo < 6 x i 0 Mpc 2 (left panel) and clusters with F500 > 6 x 10~ 5 Mpc 2 (right panel) are 
reported. The black open points are clusters with radio upper limits, while clusters with giant RH and with USSRH are shown as the red points and green asterisks, 
respectively. The vertical and horizontal dashed lines: c = 0.2 and w — 0.012. 

(A color version of this figure is available in the online journal.) 


that mergers have a crucial role in the formation of these cluster- 
wide diffuse radio sources. 

Another interesting observation is that all clusters with T500 > 
1 .3 x 1(T 4 Mpc 2 are merging clusters and host a giant RH. These 
clusters are very massive systems with M500 > 8 x 10 14 M Q . 
In particular, if we consider only clusters belonging to the 
EGRHS (plus the “Bullet” cluster lying within the same redshift 
range), we have six clusters with K500 > 1.3 x 1 0 4 Mpc 2 : 
four giant RHs and two USSRH. Why do we not find massive 
relaxed clusters in the EGRHS? The EGRHS is an X-ray 
selected sample, thus there are no reasons why we should miss 
a population of massive relaxed clusters, which are generally 
X-ray luminous. A possibility is that the F500 estimates for 
merging clusters are biased high with respect to M5 00 ■ Numerical 
simulations show that merging clusters fall below the M—Y 
scaling relation, such that their inferred masses could be biased 
low (e.g., Krause et al. 2012). However, recent observations 
based on SZ and weak-lensing cluster mass measurements show 
that merging clusters have weak-lensing masses 40% lower than 
relaxed clusters at fixed T500, so that their inferred SZ masses 
are biased high (e.g., Marrone et al. 2012). The latter authors 
suggested that the possible cause of these discrepancies could 
be found in the over-simplicity of the adopted models to fit the 
weak-lensing data. 

A more promising hypothesis is that the lack of massive 
relaxed systems in the EGRHS is due to the redshift range of this 
sample, 4 ~ 0.2-0. 4, which is not far from the formation epoch 
of these massive systems, M500 > 8 x 10 14 M 0 (e.g., Giocoli 
etal. 2007, 2012). In this case, the probability to observe massive 
relaxed clusters is smaller; we will investigate these points in 
more detail in a separate paper (R. Cassano et al., in preparation). 

7. SUMMARY AND CONCLUSIONS 

A number of correlations between thermal and non-thermal 
cluster properties, i.e., P\ 4 — Lx, P\a — M, and Pi 4 — Tx, 


have been reported for clusters hosting giant RHs since the last 
decade. However, due to the small statistics and to the lack of 
statistical samples of clusters observed at radio wavelengths, the 
reliability of these correlations and the effects of observational 
biases were not clear (e.g., Rudnick et al. 2006). Only recently, 
thanks to the GRHS (Venturi et al. 2007, 2008), it has been pos- 
sible to rely upon a solid sample of clusters with homogeneous 
and deep radio observations. For the first time, it was possible 
to place firm upper limits to the diffuse radio flux of clusters 
without extended diffuse radio emission at the detection level 
of the survey. These upper limits allowed for the study of the 
distribution of clusters in the Pi .4 — Lx and to discover a bi- 
modal behavior in the population of clusters: RH clusters lying 
on the Pi 4 — Lx correlation and radio-quiet clusters (Brunetti 
et al. 2007, 2009). Most important, the separation between RH 
and radio-quiet clusters has a correspondence in the dynamical 
state of clusters, with merging systems that harbor RHs and 
radio-quiet clusters that are statistically more relaxed (Cassano 
et al. 2010). The bimodality has been questioned in the light of 
the cross-correlation of the GRHS with the Planck ESZ cluster 
catalog. It was shown that while almost all RHs have been de- 
tected in SZ, only 4 out of 20 upper limits were detected (Basu 

2012) . This was interpreted as a weaker or absent bimodality 
in the radio-SZ plane. The proposed explanation for this was 
that SZ measurements allow an unbiased estimate of the cluster 
mass, whereas X-ray based cluster samples are biased towards 
the detection of bright cool core clusters, which may induce 
an apparent bimodal distribution of clusters in the radio-X-ray 
plane (Basu 2012). 

In this paper, we revise the radio-X-ray and radio-SZ 
correlations. Our analysis is based on the EGRHS (Kale et al. 

2013) . We searched and found information in the ROSAT and 
Chandra archive for a subsample of 40 clusters: 29 with upper 
limits to the radio powers and 1 1 with giant RHs. In addition 
to this sample, we also found information for a sample of 14 
clusters hosting well-known RHs from the literature. These are 
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used to obtain better leverage in radio/X-ray luminosities, which 
helps in the derivation of more robust scaling relations. 

First, we derive the correlation between the monochromatic 
radio power of halos at 1.4 GHz and the 0. 1-2.4 keV band 
X-ray luminosity of the parent cluster. We revaluate in a ho- 
mogeneous way the radio flux of all the halos by using GMRT 
and literature data and measure the X-ray luminosity of the 
clusters within R500 from pointed ROSAT observations and 
Chandra when ROSAT data are not available (or not sensi- 
tive enough). For the first time, we show the presence of 
a scaling Pi. 4 ghz oc /^ 0 l 0 ±0 ' 2 . Being steeper than other ha- 
los, USSRH are in general under-luminous with respect to 
this correlation. Their inclusion in the fit procedure produces 
a slightly lower normalization and an increase of the scatter 
with respect to the best-fit relation. We also correct the X-ray 
luminosity of the parent cluster by modeling the X-ray bright- 
ness distribution and excising the cool core. We find that for L500 
(or T 500, cor) ^ 5 x 10 44 erg s _1 , the distribution of clusters in 
the {Pi a ghz 1 L500) and (Pi. 4 gh z , Tsoo.cor) planes is bimodal: RH 
clusters lie on the correlation, while clusters with upper limits 
to the radio power are below the 95% confidence region of the 
best-fit correlation. This allows us to conclude that the presence 
of cool-core clusters does not affect the bimodal behavior of 
clusters in the radio power X-ray luminosity plane. 

To investigate the behavior of clusters in the radio-SZ 
diagram, we cross-checked the sample of clusters selected 
from the EGRHS with the 15.5 month PSZ catalog (Planck 
Collaboration et al. 2013b) and found SZ information for all 1 1 
RHs and for 19 out of 29 clusters with upper limits. Also, for 
the remaining 14 clusters with giant RHs, we found information 
in the PSZ catalog. We found a clear correlation between the 
RH Pi a and the cluster T500 of the form Pi .4 oc T 5 2 0 q 5±0 ' 28 , 
in line with previous findings (Basu 2012). However, contrary 
to previous findings, at least for T500 > 6 x 10~ 5 Mpc 2 
(roughly corresponding to M500 > 5.5 x 10 14 M 0 ), we find 
that all clusters with radio upper limits lie below the 95% 
confidence region of the best-fit correlation, highlighting a 
bimodal behavior of clusters in the radio-SZ diagram. This 
segregation is strengthened by the separation of those clusters 
in the morphological diagrams: clusters with diffuse radio 
emission are merging clusters, while the great majority of 
clusters with upper limits are relaxed, thus highlighting the 
importance of merging events in the generation of giant RHs. We 
also use the tight correlation between the core-excised cluster 
X-ray luminosity Psoo.nc and F500 to derive the predicted value 
of T500 for those clusters in our sample that are actually not 
detected by Planck. As expected, we found that the majority of 
them (8 out of 10) are in clusters with T500 < 6 x 10~ 5 Mpc 2 , 
where the completeness of the PSZ catalog is poor (about 20%). 
Interestingly, half of the non-detected clusters are cool core 
clusters; only two of seven cool core clusters of our sample 
were detected by Planck, suggesting that in the region of lower 
completeness, Planck loses preferentially cool-core clusters 
with respect to merging systems. 

The EGRHS is not selected in mass but in X-ray luminosity. 
However, considering that the completeness of the PSZ catalog 
for M500 ) 6 x 10 14 M 0 at 0.2 ^ z ^ 0.33 is ~80% (Planck 
Collaboration et al. 2013b) and cross-correlating the PSZ 
catalog with the EGRHS, we estimated that the completeness 
in mass of the EGRHS is ~55%, and the addition of radio 
observations of ~17 galaxy clusters from the PSZ catalog will 
provide a sample of mass selected clusters with deep radio data 
and a completeness of ~80%. For a comparison, assuming 
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the same masses and redshift range, we estimated that the 
completeness of the ESZ Planck catalog is of the order of 35%. 

Remarkably, we found that for T500 > 1.3 x 10 4 Mpc 2 
(or M500 > 8 x 10 14 M q ), all clusters of the EGRHS are in 
the process of merging and have a RH. We consider several 
possibilities to explain this result and conclude that the most 
likely explanation is that we are looking at these massive systems 
near their formation epoch (we selected clusters at z ~ 0.2-0. 4), 
and thus the probability to observe massive relaxed systems at 
these redshift should be relatively low. 

In Section 4, we derive basic scaling relations predicted by the 
two main scenarios put forward so far to explain giant RH. Under 
the assumption that synchrotron emission dominates energy 
losses of relativistic electrons in the ICM, and that the ratio 
between the energy density of cosmic -ray protons and thermal 
ICM in the radio emitting volume does not depend on cluster 
mass, “secondary models” predict that the synchrotron power of 
the halos scales as vP syn oc i (e.g., Kushnir et al. 2009) 
and vP syn oc F 5 oq 5-1 ' 43 - These scalings are flatter than those 
derived from observations in the present paper (see Table 3). 
Re-acceleration models typically predict steeper slopes. For 
example, under similar assumptions for magnetic field and 
cosmic rays, following Cassano et al. (2007), the scalings of 
the halo radio power with the cluster mass and SZ flux are 
vP syn oc Mj| 00 and vP syn oc F 5 2 qq, respectively, and are in 
agreement with the observed scalings (see Table 3). A detailed 
comparison between model expectations and observed scalings 
which considers the full range of model parameters is beyond 
the aim of this paper. 

It is also worth mentioning that in both the radio-X-ray 
and radio-SZ diagrams, clusters with USSRH are all below 
the 95% confidence region of the best-fit correlations. They are 
preferentially located in the region between “classical” RHs and 
radio upper limits. This is not surprising, since these RHs are 
steeper than those on the correlations, and thus their synchrotron 
emissivity at 1.4 GHz is lower with respect to that of RHs 
with flatter spectra. Interestingly, their position relative to the 
correlations was already predicted by models in which RHs 
are generated as a result of the turbulent re-acceleration of 
relativistic electrons in the ICM (e.g., Cassano 2010; Donnert 
et al. 2013). 
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